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I. INTRODUCTION 

Hydrodynamical models arc widely used in the descrip- 
tion of collective properties in high-energy collisions: in 
a thermal picture, the observed hadronic final state is 
treated as a result of particle production from a thermal 
ensemble, and hydrodynamics is a powerful tool to dy- 
namically connect this final state to the initial state of a 
high energy collision. The solutions of the hydrodynam- 
ical equations are thus important because they can shed 
light on the dynamics of the strongly interacting mat- 
ter. Exact solutions 1 are important also theoretically: 
they are solutions of a highly coupled set of nonlinear 
differential equations. They can be used to test numeri- 
cal codes, and in some cases, they by themselves provide 
insight into the dynamics of the hot and dense matter. 

The application of hydrodynamics in the field of high- 
energy physics has a long history. In the '50s, Lan- 
dau and his collaborators formulated relativistic hydro- 
dynamics, and gave a paradigmatic example of exact so- 
lutions, nowadays known as the Landau-Khalatnikov so- 
lution [1-3]. Another historic exact solution is the one- 
dimensional, boost-invariant accelerationless one, known 
as the Hwa-Bjorkcn solution [4, 5]. In the past years, 
the heavy-ion experiments at the RHIC particle acceler- 
ator yielded many surprising results on the the soft kine- 
matic domain of particle production in heavy ion colli- 
sions. The scaling behavior in the single-particle spectra, 
the two-particle correlations and anisotropy (flow) prop- 
erties could be explained under the assumption that the 
created matter is a strongly interacting liquid. Because of 
this, the interest in rclativistic hydrodynamics flourished 
again in the past years. 

Because of their analytic simplicity, exact solutions of 
relativistic hydrodynamics are very useful in the study 
of the space-time picture of high-energy reactions. The 



Landau-Khalatnikov solution was used for decades in 
the description of cosmic ray elementary particle colli- 
sions. The Hwa-Bjorken solution yields an estimate of 
the initial energy density of high-energy collisions [5]. 
(Ref. [6] presents an improvement of this estimation.) 
Recently, an interpolating solution between the Landau- 
Khalatnikov solution and the Hwa-Bjorkcn solution was 
discovered [7].) The Buda-Lund model [8, 9], which gives 
a good description of hadronic observables, e.g. univer- 
sal scaling of elliptic flow, relies on a class of exact so- 
lutions [10]. There arc many other exact solutions with 
important applications, see e.g. Refs. [10-12] for some 
non-relativistic, Refs. [6, 7, 13-18] for some relativistic 
examples. 

In this paper the equations of rclativistic hydrodynam- 
ics are approached from an unusual direction. There is 
a non-relativistic exact solution [12], described in detail 
in Section II, that has a remarkable property: the cor- 
responding phase-space distribution of the flowing par- 
ticles satisfies the collisionless Boltzmann equation [12]. 
This means that one can prepare a microscopic initial 
condition of the flowing particles, where the collisionless 
free motion of them maintains local thcrmalization, and 
thus the macroscopic quantities (temperature, pressure, 
velocity) solve the hydrodynamical equations (since they 
follow from kinetic theory). I will refer to this peculiarity 
as collisionless flow or Knudsen flow, in analogy to the 
case of Knudsen gases, where the collisions of particles 
play no role in the thermalization. (A hydrodynamical 
solution whose phase-space distribution satisfies the colli- 
sionless Boltzmann equation is not necessarily free of col- 
lisions, but it can be, for appropriate microscopic initial 
conditions.) The existence of such solutions may be sur- 
prising, because hydrodynamics relies on local thermal 
equilibrium, which is usually maintained by collisions, so 
it is interesting that in some cases the particles behave 
thermodynamically equilibrated even without collisions 2 . 



We make the distinction between numerical solutions (whose ad- 
vantage is generality and applicability to many different initial 
conditions) and exact ones, formulated in terms of explicit for- 
mulas, the advantage of them being the paramctrization of a set 
of initial conditions and simplicity. 



2 The mentioned non-relativistic solution [12] was discovered as 
a phase-space distribution whose general form is pertained with 
collisionless evolution, then it was realized that this is a solution 
to the hydrodynamical equations [19]. 
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In the following I investigate this idea and general- 
ize it to relativistic hydrodynamics. The result is a new 
family of exact solutions, even ones with nonzero rota- 
tion. (As far as I know, relativistic hydrodynamical so- 
lutions with nonzero curl were not known before.) These 
and other space-time characteristics of the new solutions 
make them a good candidate for using them in the de- 
scription of high-energy reactions. 

II. NOTATION AND BASIC EQUATIONS 

Let's denote the space-time coordinate by 2+ = (i,r), 
and the metric tensor by g^" = diag (1, — 1, . . .). We 
investigate 1 + 1 and 1 + 3 dimensional flows as well, 
the notation D will stand for the dimensionality of the 
space: g£ = D + 1, 8kk = D. (Greek letters de- 
note Lorentz indices, Latin letters denote three-vector 
indices.) The four- velocity of the fluid is = 7(l,v), 
where 7 = \/l — v 2 , and the thermodynamical quantities 
are: T the temperature, e the energy density, p the pres- 
sure. (Their dependence on 2+ is usually suppressed in 
the notation.) The Equation of State (EoS) of the matter 
connects £ with p and T. If we consider a fluid consist- 
ing of individual particles, then their (conserved) number 
density is denoted by n, and the corresponding chemical 
potential by u. The fundamental thermodynamical rela- 
tions are 

e + p = Ts + fin, dp = sdT + ndp,. (1) 

The hydrodynamical equations are well known, I recite 
them in the form used here. In the non-relativistic case, 
one has to introduce the particle mass tuq. The Eulcr 
and the energy conservation equations are then 

nm (5 + (vV)v)=-Vp, (2) 



_ + V( £ v) = -p(Vv), (3) 

and if a (conserved) particle number density is also con- 
sidered, we have 

^ +V (nv) = 0. (4) 

The hydrodynamical equations need to be supplemented 
by an appropriate EoS to close the set of equations. 
Throughout this paper, the following equation of state 
is used: 

£ = Kp, p = nT. (5) 

This EoS is used in most cases of exact solutions, eg. the 
Hwa-Bjorken solution, the Landau-Khalatnikov solution, 
and many others. The constant k is usually left for free 
choice, but in this paper, all the solutions investigated 



are valid only if we constrain the value of k: we must 
have tz == D/1 in the nonrelativistic case, and k = D 
in the relativistic case 3 . However, for clarity, n will be 
sometimes retained in the notation. 

In the more interesting relativistic case, the Euler 
equation, and the energy and particle number conser- 
vation equations are 

(e + p) xfdvv? = (gW - u"u p ) d pP , (6) 
(e + p)d p u p + u p d p e = 0, (7) 

9 M (nu^) = 0. (8) 

If we do not investigate the n density (and its continuity 
equation (8)), then the equations contain only i/+, p and 
£, so the e — np EoS closes the set of the equations, 
without any reference to T or n. 

Nevertheless, a central topic of this paper, the rela- 
tion between the collisionlcss phase-space evolution of the 
matter and the hydrodynamical flow, is meaningful only 
if we include a non-vanishing, conserved n, and thus /i, 
and in the spirit of Eq. (5), also T. The interesting quan- 
tity in kinetic theory, the phase-space distribution / is a 
function of the momentum p^ = (E, p) and 2+ . In the 
non-relativistic case, it has the Maxwcll-Boltzmann form 
as 

fU \ ( n (p-m v) 2 \ 

/ft,r, P ) = ex P (_- ^ T y ( 9 ) 

while the relativistic generalization is 

/( a> p) = e«p(£-*e^). (10) 

For non-vanishing n, / is required to be normalized as 

J d D pf(x,p)=n(x). (11) 

(For example, in the non-relativistic case this is fulfilled 
if exp (n/T) = n (2TrmoT)~ D ^ 2 , this is one of the reasons 
for choosing Eq. (5) to be the EoS, with the restrictions 
on k explained there.) 

There is a known class of ellipsoidally symmetric non- 
relativistic exact solutions [10]. A spherically symmetric 
special case (mentioned in the introduction) , with homo- 
geneous temperature and Gaussian density profile was 
discovered earlier [12]: 

a(t) ( a \ 3 / r 2 \ 



It must be noted that any hydrodynamical solution with e = ftp 
EoS has another important property: one can introduce the bag 
constant B, and set p— s-p — B, e— > e + B and the solutions still 
remain valid with the „bag" equation of state, which is e — B = 
K(p + B). 
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T = Tn 



a(t) 



a(t) = A(t-t y + B, (12) 



with ao, to, A, and -B constants. This solution has 
the surprising property explained in the introduction: 
the collisionless phase-space evolution [12]. (With other 
words, Ref. [12] describes a Knudsen flow.) The Boltz- 
mann equation without collisions takes a particularly 
simple form: 



9 ft \ , P 9 ti \ 
— /(r,p + — 7-/r,p 

at mo or 



0. 



(13) 



Indeed, this has the general solution / (t, r, p) = 
/ (to, r — pt/nio, p), which is a free streaming. The solu- 
tion given by Eq. (12) satisfies this, because Gaussian dis- 
tributions are stable, and the velocity field is self-similar. 
Appendix A proves that there are no essentially different 
non-relativistic collisionless solutions. 

In the relativistic case, the collisionless Boltzmann- 
cquation is 



P"0„/ = O (Vp) 



(14) 



This is the relativistic analogue of Eq. (13). In the next 
section the solutions of the coupled (relativistic) hydro- 
dynamical and collisionless Boltzmann equations are pre- 
sented. The derivations are rather technical, they are left 
for Appendix B and C. 



III. NEW EXACT RELATIVISTIC SOLUTIONS 
AND KNUDSEN FLOWS 

The equations (B14) and (B16) in Appendix B allow 
us to derive the collisionless hydrodynamical solutions 
(the Knudsen flows) . Some of them are special cases of 
known solutions, (here the new result is the collisionless- 
ness), some of them are new results in themselves. Some 
of these new solutions are special cases of new more gen- 
eral solutions (which are not Knudsen flows) . It must be 
noted that the new solutions are in a sense independent of 
the investigation presented in the previous section: they 
are valid, exact, explicit hydrodynamical solutions (p, e 
and it M fields) for the e = up (k = D) equation of state 
(as in Eq. (5)). 

The 1+1 dimensional solution given by Eq. (B16) (with 
x = x 1 being the only spatial coordinate) can be ex- 
pressed with two arbitrary functions, x( x+ ) an d £ (x~) 
(with the notation x^ = t ± x) as 



X (x+)-Hx~) 

x ( x +) + t(x-y 



T = 



To/2 



Vx(*+K(z-) 



(15) 



Recall that we must have n = D in the EoS, that is, 
here K = D = 1. In Appendix D of Ref. [18], a general 
solution of the hydrodynamical problem for k = D = 1 
was presented, with similar characteristics: the velocity 
field v and the pressure p is expressed as the combinations 



of two arbitrary wave-shapes propagating in the opposite 
direction. (The notation used there is different, the F 
and G functions introduced there are F = ln\ an( l G = 
— hi£, respectively). Eq. (15) is a special case of that 
solution, when p = nT and n = const ■ T, that is, T = 
const ■ y/p, which is a minor restriction on the general 
solution. We see now that this restriction is enough to 
have a solution with collisionless phase-space evolution. 

We now turn to the three-dimensional case. In the 
following we denote |r| by r, and make use of the Rindler- 
like coordinates r and n (,, proper-time" and ,, space-time 
pseudorapidity" ) , defined here as 



\ft? 



1, t + r 
— In . 

2 t-r 



(16) 



(Note that this notational convention favors spherical 
symmetry.) Eq. (B14) in Appendix B offers several pos- 
sibilities for relativistic hydrodynamical solutions. The 
investigation of Eq. (B14) is left to Appendix C, here 
the results are collected. The notations To, po and To 
always stand for ,, initial" values (to set the scale of the 
temperature and the pressure). 

The first case is the well-known Hwa-Bjorkcn solution 
in 1 + 1 dimensions, and the Hubblc-like (or Buda-Lund 
type) solutions [14] in D ^ 1: 



r to 
v=-, T = T — , 

t T 



P 



T 



K + l 



(17) 



The solution itself is a known result, the new result is 
that it corresponds to a Knudsen flow. 

The second case of solutions is an accelerating solution, 
which is a generalization of the previously known solution 
in Rcfs. [6, 18]. Its form is 



2tr + a 



t 2 



T 

>\To 



K+l 



T 



Tori 



(t 2 + p) + Apr 2 sinh 2 r/ — 4 (ar) t — a 2 



(18) 



with arbitrary a constant three- vector and p > constant 
(the p = 0, a = case was presented in [6, 18]). For a^O 
this new solution prefers a direction in space. The a = 0, 
p > solution is sperically symmetric, and it is finite in 
n, what is an improvement on Refs. [6, 18]. (The term 
,, finite" means here that p drops to zero at any given 
t = n hypersurface as rj goes to ±oo. Roughly speaking, 
the domain in rj where the pressure and matter density 
is different from is finite.) 

The fluid trajectories of the (already known) p = so- 
lution are uniformly accelerating in their local rest frame, 
Fig. 1 shows these trajectories. Inside the lightcone they 
describe an explosion from an initially infinitely dense 
point at rest. The trajectories of the p > solutions 
are qualitatively similar, (shown on Fig. 2 for p = 1), 
but they describe an explosion where the initial point 
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FIG. 1: Fluid trajectories of the already known [18] exact 
solution, the p = case of Eq. (18), inside (solid) and outside 
the lightcone (dashed). 

Fluid trajectories~| 
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FIG. 2: Fluid trajectories of the new, accelerating 77-fhiitc 
exact solution of Eq. (18), for the p = 1 case as illustration. 



is ,, smeared" on a distance scale of ~ ^fij. The a / 
case describes an explosion where this initial sphere has 
nonzero velocity and pressure gradient. 

We can have solutions where the initial fireball has 
nonzero rotation, so has the whole velocity field. Again 
with arbitrary constant p and constant three-vector a, 
and a new constant three- vector B, the solution is 



2ir + a + Bxr 



t 2 



P 



T 



K+l 



T 



(r 2 + p) + Apr 2 — (B x r + a) — At (ar) 



(19) 



(20) 



There exists a rotating generalization of the Hubble- 
likc how of Eq. (17) as well (for the derivation, see Ap- 
pendix C): 



r + Ef + Bxr 
t + (Er) ' 



T 



K + l 



(21) 



T = 



T to 



t 2 + (Br) 2 - B 2 r 2 - EH 2 + (Er) 



(22) 



with constant E and B three- vectors. (In a special case, 
if EB = 0, one of the E, B vectors can be set to zero.) 
This flow also describes expansion, where the perturba- 
tion with respect to the Hubblc-likc flow do not necessar- 
ily tend to zero asymptotically. (However, it is not only 
the Hubble flow seen from a rotating reference frame.) 

So far all the solutions presented are Knudsen flows, in 
the sense used in the introduction. However, if one does 
not consider the conserved particle density n, only T, p 
and one has valid new solutions (and the already 
known Hubble-like solution) to the relativistic perfect 
fluid hydrodynamical equations (not restricted to colli- 
sionless evolution of conserved massless particles). We 
see that one can have solutions describing an expansion 
of an initial fireball that can have finite size, initial ve- 
locity and directional pressure gradient as well as initial 
rotation. (To my best knowledge, rotating relativistic hy- 
drodynamical solutions were never found before.) These 
properties make these solutions themselves an interesting 
tool in the description of high energy heavy ion reactions, 
where a rotating fireball of finite size is initially present. 

If we abandon the requirement of collisionlcssness, we 
can generalize the solutions presented above. For any ve- 
locity field there exists a scalar function S which satisfies 
^| + vV5 = 0. For such S, one can take any V (S) func- 
tion, and replace T -> T/V (5), n -> n-V (S), p -> p. The 
obtained formulas solve the relativistic hydrodynamical 
equations (the only difference would be in the particle 
density continuity equation (8); it can be checked that 
the solution remains valid), but the collisionless Boltz- 
mann equation is solved only if V (S) = 1. For example, 
the S function for the solution (18) is easily found, and 
thus the following solution is a finite hydrodynamical so- 
lution with an arbitrary scaling function: 



2ir 



t 2 



S = 



To[^ 

Po 



v(sy 



no 



Pa J 



Po T o 



2(k+1) 



((' 



(23) 



- p) 2 + 4pr 2 sinh 2 n 



and since p is finite in 77 (as explained after Eq. (18)), n 
and T can be finite simultaneously. 



IV. DISCUSSION AND SUMMARY 

In this paper, the connection between hydrodynami- 
cal evolution and collisionless kinetic phase-space evolu- 
tion was explored. It turns out that in some cases the 



5 



hydrodynamical flow solves the collisionlcss Boltzmann 
equation as well: thermalization can be maintained by 
free streaming. This means that in these cases, one can 
prepare a special initial condition for the phase-space dis- 
tribution where the collisionless motion of the particles 
maintains local thermalization. 

The very early phase of a heavy-ion or hadron collision 
process, with quarks as degrees of freedom, is a domain 
where this ,,Knudsen flow" may have some implications. 
With increasing energy the lack of collisions (the decrease 
of the cross section due to asymptotic freedom) as well as 
the ultrarelativistic limit (which, from a statistical physi- 
cal point of view, is equivalent to the masslcss limit) may 
become realistic, it is interesting to see that a collective, 
hydrodynamical motion (even with acceleration) can re- 
sult from only the initial conditions for the phase-space 
distibution, without the need for collisions to achieve the 
hydrodynamical character. Thus it seems that in some 
cases the thermalization process is less connected with 
the actual collisions than one would think. 

The new exact solutions of relativistic perfect fluid hy- 
drodynamics, which emerged during the simultaneous in- 
vestigation of the collisionless Boltzmann equation and 
the hydrodynamical equations, are important on their 
own right, and are not necessarily restricted to collision- 
less evolution of conserved masslcss particles (only the 
e = up equation of state is necessary). Among these so- 
lutions, one has accelerating solutions beyond spherical 
symmetry, and realistic, accelerating solutions describing 
expansions with finite space-time rapidity profile. For 
the first time, rotating relativistic hydrodynamical solu- 
tions were also found. All these new solutions definitely 
deserve future work and investigations to establish their 
connection to the experimentally observable quantities 4 . 
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Appendix A: Nonrelativistic collisionlessness 

In this Appendix it is shown that the only reason- 
able three-dimesional non-relativistic collisionless hydro- 
dynamical solution is the one discovered in Ref. [12], and 
recited in Section II. The derivation goes in the same vein 
as in the relativistic case in Appendix B: we expand the 



collisionless Boltzmann-cquation (13), and require that 
the terms of different order in p vanish separately, for all 
p. From Eqs. (9) and (13), keeping in mind the normal- 
ization of Eq. (9), we have: 



1 / dn p _ 

- or + — V" 

n \ at rriQ 



_3_ far 

2T \~dt 



pVT\ 

m J 



(p - m v) 
2m T 2 



(p - mpv) 
T 



at m / 



dw 
~dt 



m 



We have from the terms of order p 3 and p 2 
1 dT 

VT = 0, --^S ik +diVk + d k Vi = 0, (Al) 
which has the unique solution 



T = T( 



a 2 (0) 
« 2 (*)' 



£(i) 

a(t) 



D x r, 



(A2) 



with arbitrary a (t) function and C, D constant vectors. 
Taking the curl of the Euler equation, Eq. (2), and taking 
Eq. (Al) into account, one gets that if d (t) ^ 0, then D = 
0. (The other case is not really interesting.) Eq. (A2) 
means a solution of the energy continuity equation (3) if 
and only if k = D/2. For the moment we assume C = 0, 
in that case the velocity field is spherically symmetric, 
and this implies that p, and hence n must be spherically 
symmetric. With some calculation one can now write 
down the general solution of the continuity equation for 
n, Eq. (4). What we have now is 



T = T 



o 2 (0) 
a 2 (f) ' 



a 3 (t) 



a(t) ' 



S(t,r) = 



a 2 (t) 



(A3) 



(A4) 



The v (S) function can be arbitrary at this stage. The 
variable S has the property that 5j| + vVS 1 = 0. Now we 
only have to put these in the Euler equation to obtain 
the unique solution as in Eq. (12). It indeed solves the 
Boltzmann equation, too. With some calculation one can 
be convinced that the free C vector is nothing but the 
freedom of a Galilei-transformation: we can observe this 
flow from a moving reference frame. 



Appendix B: Relativistic collisionless evolution 



4 During the finalization of this manuscript I became aware of a 
different class of new exact solution of the relativistic hydrody- 
namical equations that utilizes conformal symmetry in the trans- 
verse plane to generate radial flow [20]. This solution deserves 
future investigation as well. 



The relativistic phase-space distribution / (x,p) is 
given by Eq. (10), for a hydrodynamical flow. To have a 
collisionless hydrodynamical flow, besides the equations 
of hydrodynamics, Eq. (14) also must hold, for all allowed 

values. So, from Eqs. (14) and (10), that the first and 
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the second order terms in p 9 vanish separately, we get 
two equations. The first one is 

9^ = 0, (Bl) 

which means that /i = const. ■ T. This, together with Eq. 
(1), yields the Euler equation in a simpler form: 

IVflU" = (g" p - u^u") d P T, (B2) 

containing only the temperature and the four-velocity. 
Also for the (7) energy conservation equation, with Eq. 
(1) and the EoS from Eq. (5), we get the alternative form 
(which is essentially the expression of entropy conserva- 
tion for the case of the investigated EoS): 

0„ (TV) = 0. (B3) 

The other constraint from the Boltzmann equation de- 
pends on whether we consider a fluid of massive or mass- 
less particles. For the massive case we have 

Mf) Wt) =0 ' (B4) 

while, if they are massless (or, in other words, we consider 
ultrarelativistic particles), the constraint is weaker: 

with a free a (x) scalar function. (This weakening is 
due to the fact that for massless particles p^Pp = 0.) 
From now on, we investigate the massless case only. Con- 
tracting Eq. (B5) with u", and using Eq. (B2), we find 
Ta = —2u p d p T, so we have 

dp (£) + ft, (^) + ^u»d p T ■ 9 p V = 0. (B6) 

This equation contains the Euler-equation (we get it by 
contracting with u v ), and contracting Eq. (B6) with g^ v 
yields dp (T ^) = 0, which is essentially the entropy 
conservation equation, since we have n = D. (This is 
a reason for using ft = D in the relativistic case.) It 
is clear from the text before Eq. (B3) that this form of 
the entropy conservation equation is equivalent to the 
energy conservation law, Eq. (7), so we see that Eq. (B6) 
is the only equation we have to solve. For this purpose 
we introduce the four-vector = (<p,A.) = u^/T. The 
three- velocity is v = A/cp. With this, Eq. (B6) becomes 

APA° {2gp U dpA a - g p<J (3pA v + 8„Ap)} = 0. (B7) 

With some investigation (e.g. writing this equation in 
non-relativistic notation) we find that Eq. (B7) is equiv- 
alent to the one we get by dropping A 9 A a : 

9p V (d P A a + d a A p ) - g pa (dpA v + d v Ap) = 0. (B8) 

Remarkably, we have obtained a set of linear differential 
equations. For solving this, let's expand Ap in a Taylor- 
like series: 

Ap = a^+a$x»+al%x»x 9 +a$ p!T x» X Px° + . . . . (B9) 



The coefficient tensors are by definition totally sym- 
metric in their last i — 1 indices. This expansion works 
if both T and vf* arc smooth, and T / 0. (Wc should 
have such domain in space-time if we are after hydrody- 
namical solutions.) From the requirement that Eq. (B8) 
holds order by order in x M , we have 

(2) _ (2) (3) _ (3) 

9p<r a [pv] — 5V a [p CT ]' 9p<J a [p U \\ — 5^ a [ pfT ]A> 

9p^l ]Xv = g^a\pl ]Xv , (B10) 

where square brackets mean symmctrization. To pro- 
ceed, we must now distinguish between the D = 1 and 
the D 1 case. 

In the D ^ 1 case, Eq. (B10) is satisfied for ajxv if 

a pv = 19p.v "F Fftv , 

with arbitrary given 7 constant, fo M vector and F^" anti- 
symmetric tensor. The equations for the a^s for i > 2 
imply 

«^...+«2a...=^, CBll) 

with totally symmetric 7a... tensors. Taking the antisym- 
metric part of this equation in fi and A, and then in fx 
and p, using the symmetry properties of the a^s, and 
again Eq. (Bll), we obtain 

a?L. = \ (V^f. + 9,xT^. - g vX T$) . (B12) 

(3) 

For i = 3, taking — —6a, wc immediately have the 
result as 

%v\ = 2 (9vx b » + 9i*vb\ ~ 9u\bp) (B13) 

For i > 3, that o'^ be symmetric in A <H> 77, we can 
express these constraints as 

(v^t.-^...) - (W^ll. -ft,^...) = 0. 

For every set of the other remaining indices (denoted 
collectively by dots) this homogeneous linear system of 
equations has D (D + 1) (D (D + I) + 2) /8 independent 
components (because the antisymmetry in \x f-> A and 
v <-> 77, and the symmetry in nX <-> 1/77), while a symmet- 
ric tensor T M „ has (D + 1) (D + 2) /2 independent com- 
ponents. So for D = 2 there are 6 equations for 6 un- 
knowns, for D > 3 there are more equations than un- 
knowns. Furthermore, this set of linear equations is non- 
singular, so there can be no nontrivial solutions for i > 4 
and D > 2. Combining these together, we arrive at the 
general solution for Eq. (B8) as 

t T v \ up 

A 9 = a v+ 1 x f * + F f * u x„ + (b u x„)x> 1 -± ^— , (B14) 
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where dp = a}, is an arbitary four-vector. From this, we 
get the solutions derived in Appendix C and described in 
Section III. 

The D = 1 case case is simpler: Eq. (B8) yields only 
the following conditions for ip (= A°) and A = A : 

OA = dip OA = dp 
dt dx ' dx dt 

(Here x = x 1 is the only spatial coordinate.) Again with 
the notation x^ = t ± x, this has the following general 
solution, with arbitrary \ ( x+ ) an d £ (%~) functions: 

A= = , ¥= = , (B16) 

Jo 

where an arbitrary temperature scale To is introduced. 
Having the definition of A^ in mind, we arrive at the 
solution as in Eq. (15). 

We thus found all those relativistic hydrodynamics! so- 
lutions, where thcrmalization can be maintained without 
the collision of particles (in the case of massless parti- 
cles). 

Appendix C: Three-dimensional collisionless flows 

This Appendix investigates the solutions stemming 
from Eq. (B14). We use the Rindler-like coordinates in- 
troduced in Eq. (16). 

• If the vector and the F^ v tensor in Eq. (B14) 
is zero, then, if 7 = 0, the solution is a medium at 
rest, but if 7 7^ 0, then et M can be eliminated with 
a translation x M — > x 11 — 0^/7, and introducing the 



notation To To = I/7 (to set the scale of the initial 
temperature, we obtain the solution cited in Eq. 
(17), the Hubble- like expanding solution. 

• We can have 6 M 7^ in Eq. (B14), and we can re- 
define F^ v and a M in Eq. (B14) and make a trans- 
lation in x^ to make 7 vanish. The simplest case 
then is when 6 M is time-like and F^ iV = 0: we can 
adopt a frame where 6 M = (C,0), and we obtain a 
new accelerating solution (with a M = (C/f/2, C a /2) 
and 2/C = T r 2 ) as in Eq. (18). 

• For time-like 6 M and nonvanishing we obtain 
accelerating flows with curl. We can choose a frame 
where 6 M = (C, 0), and adopt the three-dimensional 
notation for the F^ tensor as Foi = hC^i, Fik = 
\C e ikiBi, take a M = (Cp/2,£a/2). It turns out that 
with an appropriate translation and redefinition of 
the other constants, we can eliminate E. We thus 
have the solution as in Eqs. (19) and (20), with 
7W = 2/C. 

• The rotating generalizations of the Hubblc-flow, 
given in Eqs. (21) and (22) can be derived from 
Eq. (B14) if we take f = 0, 7 / 0, so we have 
A p. = Jxt i + p^ Xv ^ with the notation F oi = jE u 
Fik = T&iklBi, the solution is given by Eqs. (21) 
and (22) (with T r 7 = 1). 

• For space-like 6 M , we can choose a reference frame 
where = (0, b), and obtain collisionless ,, solu- 
tions" . However, they are not physical, since the 
velocity of some initially slow fluid trajectories will 
approach and surpass the speed of light. 
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